Cavity-QED models of switches for attojoule-scale nanophotonic logic 
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Quantum optical input-output models are described for a class of optical switches based on cavity 
quantum electrodynamics (cavity QED) with a single multilevel atom (or comparable bound system 
of charges) coupled simultaneously to several resonant field modes. A recent limit theorem for 
quantum stochastic differential equations is used to show that such models converge to a simple 
scattering matrix in a type of strong coupling limit that seems natural for nanophotonic systems. 
Numerical integration is used to show that the behavior of the pre-limit model approximates that of 
the simple scattering matrix in a realistic regime for the physical parameters, and that it is possible in 
the proposed cavity-QED configuration for low power optical signals to switch higher-power signals 
at attojoule energy scales. 

PACS numbers: 42.79.Ta,42.50.Pq,42.50-p 
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It has long been appreciated that in cavity quantum elec- 
trodynamics (cavity QED) with strong coupling [3, Q], 
the transmission or reflection of an optical field coupled 
to the cavity mode can be controlled by the state of 
a single intra-cavity atom (or comparable bound sys- 
tem of charges, e. g., as in solid-state systems) 0, 0, 
This simple insight leads naturally to various schemes 
for the implementation of optical switches in which a low- 
power beam is used somehow directly to manipulate the 
atomic internal state, and thus to determine whether a 
higher-power beam incident upon the cavity input cou- 
pling mirror is transmitted or reflected (for related pro- 
posals see d, 0, As cavity QED switches should be 
realizable in integrated nanophotonic platforms 0, [l(| 
and could potentially function down to the single photon 
level [ll| (although it is not clear whether robust circuit- 
level operation could be achieved with few-photon signals 
in the presence of finite propagation losses) this type of 
device could contribute substantially to the development 
of ultra-low power nanophotonic signal processing, for 
example in the context of on-chip nanophotonic inter- 
connect m m m . 

In quantitative performance studies of cavity-QED 
switch designs it is important to consider fundamental 
models that include the effects of optical shot noise, 
spontaneous emission, and dipole fluctuations. Such 
models should be cascadable, and should clearly indi- 
cate the key parameters for optimizing switching per- 
formance. The purpose of this article is to describe 
quantum-optical input-output models for an elementary 
class of cavity QED switches, which can be cascaded 
easily using the series and concatenation products for 
quantum stochastic differential equations (QSDE's) |l5T |. 
and to show that these models converge to a very simple 
scattering-matrix description in a certain strong coupling 
limit that seems natural for nanophotonic systems. The 
limit model may provide a useful abstraction for circuit- 
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FIG. 1: Detail diagrams for a set-reset flip-flop switch 
(adapted from [l(J): (a) signal connection diagram, (b) phys- 
ical diagram of optical inputs and outputs, and (c) required 
levels and couplings for the intra-cavity 'atom.' 



level analysis and synthesis of photonic signal processing 
systems 0> O EM IS HI] i while numerical integration 
of the primary model can be used to characterize devia- 
tions from ideal switching behavior as well as the limits 
of low-power operation for finite values of the physical 
parameters. 

Fig. [1] presents the basic structure of the cavity-QED 
switch we will consider, which was recently proposed 
in The device functions essentially as a set-reset 

flip-flop switch. The intra-cavity atom, whose internal 
state determines the overall state of the switch, is as- 
sumed to have levels {\g), \h), |e), |s)} in the arrangement 
shown. States {\g}, \s), |e)} are assumed to have the same 
angular momentum m (eigenvalue of angular momentum 
along the z-axis) while \h) should have angular momen- 
tum m — 1. Although we have drawn |e) and \s) at dif- 
ferent energies, we note that this would not necessarily 
have to be the case if other selection rules (associated 
with a symmetry, rather than energy) could be used to 
prevent the POWER cavity mode from coupling to the 
\g), \h) <-> \s) transitions. In normal operation the atomic 
state should be either \g) or \h), as |e) is ideally never 
populated and \s) serves only to facilitate transitions be- 
tween | g) and \h). We assume that the POWER input 
drives a cavity mode that couples only to the \g) «-> \e) 
transition, that the SET input drives a cavity mode that 
couples only to the \g) <-> \s) transition, and that the RE- 
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SET input drives a cavity mode that couples only to the 
\h) *-* \s) transition. While the flexibility of, e.g., pho- 
tonic crystal resonator design could potentially provide 
the three required atom-field couplings via three modes 
in a single nanophotonic structure, we here propose a 
configuration of two Fabry-Perot-type cavities to illus- 
trate concretely that the selection rules we have assumed 
are reasonable. Specifically, we assume that the cavities 
are sufficiently birefringent [2l| that the cavities support 
linearly-polarized resonant modes. In the POWER/SET 
cavity we utilize modes that are polarized along the pr- 
axis (atomic quantization axis), so that these fields only 
induce n transitions. In the RESET cavity we utilize a 
mode polarized along the x-axis, so that the associated 
field does not couple to n transitions but can couple the 
<r+ transition \h) <-> \s). We assume that the cavities are 
symmetric (as opposed to single-sided). 

We specify our primary switch model via the (S, L, H) 
coefficients [15( for a (right-sided) QSDE with S = /, the 
components of the vector L given by 

Li = L 2 — k ly /n^a, L 3 = L 4 = k 2y /K^b, 

Lb = L 6 = k 2 y/K^c, L 7 = Vfcr gs , L 8 — VTa hs , 

L 9 = VTdge, L 10 = VT<T he , (1) 

and 

H = iklg p (a^a ge - aal e ) +ik 2 g s (b ji crg S -bal s ) 

+ik 2 g r (c 1 <r hs - arjj), (2) 

where a, b and c are annihilation operators for the 
POWER, SET and RESET cavity modes, a ge = |fif><e|, 
(T he = \h)(e\, <7 gs = \g)(s\, and a hs = \h)(s\. The pa- 
rameter r is an atomic spontaneous emission rate, n p s r 
and g PtSl r are the field decay rates and atomic coupling 
strengths (vacuum Rabi frequencies) of the POWER, 
SET and RESET cavity modes, and ki. 2 are dimension- 
less scaling parameters that we will use to derive the 
scattering-matrix limit. Note that the rates of the four 
atomic spontaneous emission processes could vary with 
minor impact on switch performance. With the above or- 
dering of the components of L we specify QSDE modes 
1 and 2 as the POWER input and output, modes 3 and 
4 as the SET input and output, and modes 5 and 6 as 
the RESET input and output. 

The operating principle of the switch is as follows. 
With the atom in the \h) state, the POWER cavity 
mode does not couple to the atom at all and therefore 
the POWER input field is routed to the OUT output 
(transmitted through the cavity). With the atom in the 
\g) state, however, the POWER cavity mode experiences 
a large vacuum Rabi splitting and the POWER input 
is therefore routed to the OUT output (reflected from 
the cavity). When the atom is in the \g) state, injec- 
tion of photons into the SET cavity mode induces Rabi 
oscillation on the \g) <-> \s) transition that is promptly 
terminated by decay into the atomic \h) state via cavity- 
enhanced (via coupling to the RESET mode vacuum) 
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FIG. 2: Numerical integrations of the pre-limit master equa- 
tion, displaying the intra-cavity field amplitude (normalized 
by 0) as a function of time. In each subplot, the black (bro- 
ken) curve is for initial state \g) and the red (solid) curve is 
for initial state \h). See text for parameters. 

spontaneous emission. When the atom is in the \h) state, 
decay into the \g) state can likewise be induced by injec- 
tion of photons into the RESET cavity mode. We assume 
that under normal operation photons are never injected 
into both the SET and RESET cavity modes simulta- 
neously, as this would induce a 'race' condition which 
should not normally occur in a classical logic circuit. 
Note that when no photons are injected into either the 
SET or RESET cavity mode, the atomic state should ide- 
ally remain constant ('hold' condition), although again 
this should not be necessary in conventional circuits. 

In order to prove that this primary model should be- 
have according to these principles in a regime of strong 
coupling we invoke the QSDE limit theorem of [22| . 
first taking the scaling parameter ki — > oo and subse- 
quently taking k 2 — >■ oo. For the first limit we satisfy 
the structural requirements of [22j by choosing the limit 
state-space Hq = span{|g0 a 7i{,n c ), |/i a nf,n c ), |s0 a n;,n c )} 
(where the photon numbers n a .b,c indicate corresponding 
Fock states of the cavity modes), and for the second limit 
we choose Hq = span{|g0 a b c ), |/i0 a b c ), |s0 Q b c )}. 
In the resulting model which is actually still an inter- 
mediate model from which we will subsequently derive 
an even simpler scattering matrix — the |e) excited state 
and all three cavity modes have been adiabatically elimi- 
nated. As a result we can drop modes 9 and 10 from the 
QSDE model, as they were used only to provide spon- 
taneous emission from |e). The remaining coefficients of 
the 'intermediate' QSDE are 

Sn = S 22 = II S , Si 2 = S 2 i = —Hhs, 
S34 = S43 = $56 = Se5 = Sj7 = Sss = 1, 

L 7 = Vfa gs , L$ = VT<7 hs , (3) 
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with H and all remaining components of S and L equal 
to zero. Here II g = \g}{g\, Hh s = \h)(h\ + \s)(s\, and we 
have defined a new parameter 7 = g 2 / k (which we will 
later limit — ► 00) assuming g s = g r = g and k s = n r = n. 

We note that in the context of a circuit analysis the 
series product for QSDE's [HI can be used to connect the 
inputs and outputs of our switch model to those of other 
components. For illustrative purposes here we simply ap- 
ply coherent-state inputs of fixed amplitude /?, a s and a r 
to the POWER, SET and RESET inputs (note that in 
this type of QSDE model, coherent amplitudes are nor- 
malized such that their square-magnitudes correspond to 
photon flux per unit time). This results in the simple 
modifications L 1— > L + Sd and H 1— > H + ImjLtSd}, 
where d is a column vector of displacement amplitudes 
(di = f3, e?3 = a s , d§ = a r and all other components 
zero). The master equation corresponding to a QSDE 
model is in general 
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which now takes the form 



Pi 



-i[H', p t ] + (r + 2 7 ) ( o gsPt a\ s - -{IT, Pt ) 
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+(r + 2 7 ) ^a hs p t al s --{u s ,p t } 
+|/?| 2 (n gPt n g + n hsPt u hs - Pt ) , 



(5) 



with 



H ' = iVl( a s a ls - a * s Vgs) + iVli a r^hs ~ K^tis)- (6) 

Note that this H' includes terms from originating from 
the Lindblad part of Eq. (01 . It is easy to see that the 
terms in H' induce the atomic Rabi oscillations men- 
tioned in our description of the switch operating prin- 
ciples above, while the Lindblad terms provide the re- 
quired cavity-enhanced spontaneous emission (with rate 
r + 27) from \s) to \g), \h) as well as dephasing (at rate 
|/3| 2 ) of \g) relative to the subspace spanned by \h) and 
|s). The dephasing occurs since an outside observer can 
easily distinguish which of these subspaces the atomic 
state is in, simply by measuring the optical power in the 
OUT and OUT outputs, and is generally consistent with 
the desired switch dynamics. It would appear, however, 
that for \(3\ ^> \a r \ there should be some 'Zeno-like' sup- 
pression of the reset action which can limit the degree of 
signal regeneration |/3|/|a r)S | that can be accomplished 
by this type of switch without loss of switching speed. 
When a s = a r = (hold condition) we have relay equi- 
librium states \g) and \h). Looking at the upper 2x2 
block of S we find that it becomes the identity for \g), 
indicating perfect reflection of the POWER input beam 
into OUT, while for \h) it indicates transmission of the 
POWER input into OUT. With a s = 0, a r ^ (RE- 
SET condition) only \g) is an equilibrium state while with 
a s =/= 0, a r = (SET condition) only \h) is an equilibrium 



state, so we find clear agreement between the structure 
of this intermediate model and the switch operating prin- 
ciple described above. 

We now take a final limit 7 — > 00 to obtain a very 
simple scattering matrix model of an idealized cavity- 
QED switch. Returning to the un-driven intermediate 
QSDE model of Eq. §5§ and using once again the limit 
theorem of [22| with H = span{|t/0 a 0;,0 c ), |/i0 a 0fc0 c )}, 
we obtain the limit QSDE coefficients 

Sn = S22 = n„, S12 = S21 = —lift, 
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S44 = — ILy/2, S34 = S43 = 1 — Ug/2, 

S&q = — II/j/2, S56 = Se,5 = 1 — IIft/2, 

S36 = — 5*46 = —<Tgh/2, 

S54 = Sq3 = 5*64 = — Cftg/2, (7) 

where U h = \h)(h\, a hg = \h)(g\, and a gh = \g)(h\. All 
remaining components of S as well as both L and H are 
zero. If we once again consider displacements of input 
modes 1, 3 and 5 with coherent amplitudes (3, a s and a r 
respectively, the resulting master equation for p t can be 
written in terms of independent matrix elements in the 
{\g),\h)} basis, 

Pgg = \ (-\ a s\ 2 Pgg + \a r \ 2 P hh) , 

Phg = -^ s a;~ P/lg (|/3| 2 + |a s | 2 /2 + |a r | 2 /2).(8) 

We thus find that the rates of SET/RESET action are 
given in this limit by |a SiI .| 2 /2. At equilibrium we find 

Pgg -> K| 2 /(|a s | 2 + K| 2 ), 



Phg 



(9) 



which shows that in fact sensible operation can be ex- 
pected even with a s a r ^ (race condition). 

As mentioned in the introduction, it is straightforward 
to perform numerical integration of the primary (pre- 
limit) model given in Eqs. (TTJ) and © to characterize 
the switch performance for any desired set of finite val- 
ues of the physical parameters. In Fig. [5] we display the 
results of such an integration, performed using the Quan- 
tum Optics Toolbox for Matlab 23], for parameter val- 
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and a Sir as indicated in each subplot (with fixed 
at unity). Assuming that g values generally can be de- 
creased and k values increased via resonator design opti- 
mization, this hierarchy of rates is consistent with the pa- 
rameters (g, k, r)/27r = (16, 16, 0.1) GHz of current work 
with GaAs nanophotonic resonators and InAs quantum 
dots [ll| (although quantum dots would not seem to have 
the multi-level structure required for our scheme), and 
with projected numbers (g, k, T)/2-k = (2.25,0.16,0.013) 
GHz for GaP nanophotonic resonators and diamond-NV 
centers [Toj | . 

From the usual cavity QED input-output rule for two- 
sided resonators [HI, Hj] we have that the coherent ampli- 
tude of the reflected OUT mode will be 
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FIG. 3: Configuration of two elementary switches to match 
input and output wavelengths, in case this cannot be done 
at the individual switch level: (a) port connections, and (b) 
modified coupling diagram in switch R' (note that the atomic 
level structure is the same in R and R'). 
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FIG. 4: Detail diagrams for a set-reset double-pole double- 
throw (DPDT) switch: (a) optical inputs and outputs, (b) re- 
quired level diagram and couplings for the intra-cavity 'atom.' 

that of the transmitted OUT mode will be The 
results in Fig. [2] thus clearly show that when the SET in- 
put is active the relay is driven to the \h) state in which 
OUT is high, whereas when the RESET input is active it 
is driven to the \g) state in which OUT is high. In the lat- 



ter condition a finite contrast (power) ratio OUT/OUT 
« 66 results from the finite ratios 9p,s,r/(^p,s,r^), but 
assuming this is acceptable the results show that switch- 
ing is achieved with a power 'gain' of |/3| 2 /|a S!r | 2 = 10. 
We furthermore see that the switching occurs within a 
timescale r ~ 300 and is thus effected by r|a S!r | 2 <~ 10 
photons, corresponding to a switching energy in the aJ 
range assuming visible or near-IR wavelengths. 

While we have proposed a concrete model in which 
the POWER beam differs in wavelength from the 
SET/RESET beams (in order to utilize an atomic model 
with simple selection rules) , it should be possible to cre- 
ate circuits from such components in which many sequen- 
tial logic operations are performed with only two distinct 
signal wavelengths. To illustrate this point, we sketch 
in Fig. [3^ a cascade of two switches that achieves sig- 
nal power regeneration without overall change in wave- 
length. The first switch R is of the type we have con- 
sidered above while the second switch R' is an analo- 
gous device constructed with the modified coupling dia- 
gram shown in Fig. [3)3. We likewise note that our ba- 
sic design can be extended straightforwardly to design 
more complex switches. In Fig. [4] for example we illus- 
trate the basic principle of a device in which a single 
pair of SET/RESET input beams simultaneously switch 
two POWER beams with different wavelengths [16| : the 
POWERl and POWER2 beams are assumed to probe 
distinct longitudinal modes of the cavity. 
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